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Abstract. The performed systematic meta-analysis of the quality of data description (QDD) of existing 
event generators of nucleus-nucleus collisions allows us to extract a very important physical information. 

Our meta-analysis is dealing with the results of 10 event generators which describe data measured in the 
range of center of mass collision energies from 3.1 GeV to 17.3 GeV. It considers the mean deviation 
squared per number of experimental points obtained by these event generators, i.e. the QDD, as the 
results of independent met a-measurements. These generators and their QDDs are divided in two groups. 

The first group includes the generators which account for the quark-gluon plasma formation during 
nuclear collisions (QGP models), while the second group includes the generators which do not assume 
the QGP formation in such collisions (hadron gas models). Comparing the QDD of more than a hundred 
of different data sets of strange hadrons by two groups of models, we found two regions of the equal 
quality description of data which are located at the center of mass collision energies 4.4-4.87 GeV and 
10.8-12 GeV. At the collision energies below 4.4 GeV the hadron gas models describe data much better 
than the QGP one and, hence, we associate this region with hadron phase. At the collision energies 
between 5 GeV and 10.8 GeV and above 12 GeV we found that QGP models describe data essentially 
better than the hadron gas ones and, hence, these regions we associate with the quark-gluon phase. As 
a result, the collision energy regions 4.4-4.87 GeV and 10.8-12 GeV we interpret as the energies of the 
hadron-quark-gluon mixed phase formation. Based on these findings we argue that the most probable 
energy range of the QCD phase diagram (tri)critical endpoint is 12-14 GeV. The practical suggestions 
for the collision energies of the future RHIC Beam Energy Scan Program are made. 

BAGS numbers: 25.75.-q, 03.65.Pm, 03.65.Ge, 61.80.Mk 


I. INTRODUCTION 

The correct determination of the threshold energy of Quark Gluon Plasma formation (QGP) is one 
of the major goals of modern heavy ion physics. The experiments planned at FAIR [T] and NICA [2] 
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are targeted to study the properties of exotic nuclear matter created in nucleus-nucleus (A-|-A) collisions 
and to figure out a precise position of the critical endpoint of the deconfinement phase transition [3]. 
The Beam Energy Scan program performed at RHIC revealed somewhat different properties of strongly 
interacting matter created in A+A collisions at energies below and above y/sisiN = 11-5 GeV [1]. At the 
same time, the experiments on central heavy ion collisions done at AGS and SPS found such remarkable 
irregularities as the Strangeness Horn in A+/7r+ ratio at ^ysNN = 7.6 GeV [5] and the peaks of A/p and 
A/tt~ ratios at energies about ^/sNN = 5 — 6 GeV [6l[7]. Although there exist some claims [Hl-fTO] that the 
Strangeness Horn signals about the onset of deconfinement, the convincing explanation of its relation to 
the deconfinement is, in fact, absent. Despite the fact that a very successful description of the and 

A/tt” peaks is achieved recently in |111413] . such a success does not provide an apparent relation between 
the strangeness enhancement at these collision energies and the deconfinement transition. 

Moreover, nowadays not only the situation with the “signals” of QGP formation is unclear, but also 
there is no consensus on its formation threshold energy. For example, the authors of Strangeness Horn 
“signal” [H], claim that the onset of deconfinement starts at y/s^N = 7.6 GeV jini. On the other hand, 
the authors of [TTj state that the onset of deconfinement occurs at = 6.3 GeV, while very recently 

there appeared the works [la [IB] which are arguing that the hadron-QGP mixed phase is formed at 
y/sNN = 4.2 — 4.9 GeV. Regarding these threshold energy values, we would like to stress that the works 
[HHiniiiii do not provide the underling statistical or hydrodynamic or microscopic models which are able to 
connect their “signals” with the very fact of the hadron-QGP mixed phase formation. In Refs. daiBI the 
underlying hydrodynamic model of phase transition which explains the observed irregularities is developed, 
but similarly to the works mentioned above, i.e. [IS HD, the model of Refs. [laiiBi is also relying only on 
a certain set of experimental data from many sets of existing data. 

At the same time various sets of experimental data (the transverse mass or transverse momentum 
distributions, the longitudinal rapidity distributions of several hadron species, the total or/and midrapidity 
yields of various hadrons e.t.c.) are, in principle, available and the results of a few event generators of A-|-A 
collisions are published. These generators can be divided into two major groups: the first group includes 
the QGP state in their description (QGP models), while the second group does not include the QGP 
formation, i.e. Hadron Gas (HG) models. If the existing generators of A-|-A collisions contain just a grain 
of truth, then, we believe, it is worth to systematically compare their ability to reproduce experimental 
data and to find out this grain from such a comparison. Hence, in this work we suggest to employ both the 
available sets of data and the theoretical models describing them in order to elucidate the correct value of 
the threshold energy of the hadron-QGP mixed phase. 

Our main idea is based on the assumption that the HG models of A-|-A collisions should provide worse 
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description of the data above the QGP threshold energy, whereas below this threshold they should be able 
to better (or at least not worse) reproduce data compared to the QGP models. Furthermore, we assume 
that both kinds of models should provide an equal and rather good QDD at the energy of mixed phase 
production. Hence, the mixed phase production threshold should be slightly below the energy at which 
the equal QDD is changed to the essential worsening of QDD by HG models. 

The second primary aim of this work is to make a comparative study of two classes of the event generators 
of A+A collisions. This is necessary to hx the present days status quo, to make a realistic plan of further 
experiments on A+A collisions and to formulate the most important tasks for theoretical studies. 

The work is organized as follows. In the next section we present the basic elements of our meta-analysis. 
Section III is devoted to a detailed description of two ways of averaging used in this meta-analysis on the 
example of the data sets available at ^/s]^ = 4.87 GeV. The results and their interpretation are given in 
section IV, while section V contains our conclusions and practical suggestions for planning experiments. 

II. ANALYSIS OF DATA AND MODELS 

As it was discussed above, the existing approaches [iniiii! suggest that the onset of deconfinement 
begins somewhere between the highest AGS and the lowest SPS energies, i.e. in the collision energy range 
y/sNN = 4.2 GeV [121 dn] " 7.6 GeV [SHTU]. The main reason for such a range is that almost all irregularities 
observed either at kinetic freeze-out [HHIQ] or at chemical freeze-out belong to this energy range. 

Therefore, to compare the QGP and HG models we choose the AGS data measured in Au+Au collisions 
at ^ysNN =3.1 GeV as the lower bound of collision energy. On the other hand, we extended the upper 
bound of collision energy from = 7.6 GeV to y^sV/v =17.3 GeV, since there are arguments |T8] 

that at y/sNN ~ 10 GeV there may exist the tricritical endpoint of the strongly interacting matter phase 
diagram, which in |18j is mistakenly called a triple point. Hence, we will be able to investigate the wider 
energy range than the most interesting one. 

Our choice of collision energies y'+vlv = 3.1, 3.6, 4.2, 4.87, 5.4, 6.3, 7.6, 8.8, 12.3 and 17.3 GeV is also 
dictated by the fact that exactly for these energies there were multiple efforts to describe the experimental 
data both by the QGP and by the HG generators. Therefore, in the chosen collision energy range a 
comparison of the QDD provided by these two types of models can be made in details. The collision 
energies y^+vw < 4.87 GeV correspond to Au+Au reactions studied at AGS. At y^+vlv = 5.4 GeV 
the reactions Pb+Si, Si+Si and Si+Al were also studied at AGS, while higher values of collision energy 
correspond to Pb+Pb reactions investigated at SPS. 

The main object of our meta-analysis is the mean deviation squared of the quantity of the 
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model M from the data per number of the data points Ud for a given particle type h 
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where jg ^n experimental error of the experimental quantity and the summation in Eq. 

Q runs over all data point at given collision energy. To get the most complete picture of the A+A 
collision process dynamics, we have to compare the available data on the transverse mass (tti't) distributions 
^ ~ mT dn^dy^^ ’ longitudinal rapidity (y) distributions A = and the hadronic yields (Y) 

measured at midrapidity A = or/and the total one, i.e. measured within dvr solid angle, since right 

these observables are traditionally believed to be sensitive to the equation of state properties Hi- 


Based on the theory of measurements |20j . we consider the set of quantities 


the meta-measurements of the same meta-quantity, i.e. the QDD of j data, with the different meta¬ 

devices, i.e. models {M}, with the hadronic probe h. The whole set of quantities {x^ln)\ allows one 

M 

to find the mean value {x^In)M properly averaged over the experimental data and over models belonging 
to the same class, i.e. M=HG or M=QGP, and over all hadronic species. Using the averaged values 
for two classes of models {x^/n)HG and {x^/n)QGP, we could find the regions of their preferential and 
their comparable description of the experiments. The latter case could provide us with the most probable 
collision energy range of the hadron-QGP mixed phase threshold. Unfortunately, the published articles 


as the results of 


M 


usually do not provide one with the quantities . An exception is the article [21], in which one 

M 

can find the desired values for the m-r spectra of A-hyperons for 4 rapidity intervals measured in Au-|-Au 
collisions at highest AGS energy. Therefore, using the modern software we calculated the set of quantities 


ix^/n) 
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and their errors S{x^ jn)\ 


from the published papers. 


M 


The theoretical models taken for the present analysis belong to two groups: 


• The used HG models are as follows: ARG |22|; RQMD2.1(2.3) |53|, HSD j24LI25| . UrQMDl.3(2.0, 
2.1, 2.3) |26] . statistical hadronization model (SHM) |27j and AGSHIJET_N* |28l 129] . These models 
do not include the QGP formation in the process of A-|-A collisions. The results of the HG models 
were taken from the following publications: ARG |2T1[29|, RQMD2.1(2.3) [211 [30], HSD [SHIHI], 
UrQMDl.3(2.0, 2.1, 2.3) |31L 132] . SHM |32L [33] and AGSHIJET_N* |29j . Eurther details on what 
data at what energies were analyzed are presented in Tables I and HI-XL 


• The used QGP models are as follows: Quark Gombination (QuarkGomb) model |34j . 3-fluid 
dynamics (3FD) model [35l[37] . PHSD model [381E9] and Gore-Gorona model [10l[lT]. These gener¬ 
ators explicitly assume the QGP formation in A-|-A collisions. The results of the QGP models were 
taken from the following publications: QuarkGomb model |34j . 3FD model |36[ [37], PHSD model 
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and Core-Corona model 
Tables I and III-XI. 


More details on the analyzed data and energies can be found in 


A short description of these models along with the criteria of their selection can be found in the Appendix. 


Our meta-analysis requires the well-defined errors for the quantities ln)\ 
according to the rule of indirect measurements m as 
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where in deriving the second equality above we calculated the partial derivatives using Eq. 0 and then 
applied Eq. 0 once more. 

In order to thoroughly estimate a correspondence between the experimental data and their model 
description it is necessary to have very detailed experimental data which cover rather wide kinematic region 
and include many hadronic species. In practice, however, the available experimental information is rather 
limited and, additionally, its comparison with theoretical models in many cases is done not for all available 
data, but for certain sets only. Therefore, first of all we restricted our probes to the strange particles which 
include charged kaons K^, and (j) mesons, and also A(-|-S^), A, and hyperons. This choice was 
dictated by the fact that strange particles are the “clean” probes, since they are created at primary/hard 
collisions. As it was mentioned in the Introduction, several existing “signals” of deconfinement transition 
are based on the characteristics of AT-mesons [HHio], hence it was natural to consider the strange particles 
first. 


Then for a given probe h and an observable A we calculated the average of {x^/n)\ 
of the same class as 


over the models 
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where the symbol M defines the class of models, i.e. M G {HG;QGP}, which are averaged with the 
weights (jo{M). Here Nm is the number of used theoretical models. 

In order to verify the stability of our hndings we employed two ways of averaging in Q. Eirst of them 
is an arithmetic averaging with the equal weight for all models, i.e. 
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where Nm is the upper limit of the sum in Q, or the number of used models. For the arithmetic averaging 
we calculated the averaged errors according to the formula 
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For the weighted averaging Q the best estimate for the average error is given by 
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same quantity and if the corresponding errors A^ix"^/N)\ 

In our case the measured quantity is the quality description of the observable A. 

Each way of averaging has its own advantages. Thus, the arithmetic averaging with the weights Q 
includes all measurements on equal footing and, hence, as we will see, it allows one to equally account 
for the contributions coming from different kinematic regions. On the other hand, the weighted averaging 
of Eq. Q ‘prefers’ the measurements with the smallest value of the QDD error and it provides the best 
estimate for the measured quantity [50], i.e. for QDD in our case. Note that we also used alternative 
ways to average the QDD, but usually found the results similar to one of two ways of averaging used 
here. Therefore, we concentrate on the averaging methods given by Eqs. and Q, since they are more 
convenient than the other ones and they have a well-defined meaning within the theory of measurements 


These two ways of averaging are used further on for averaging over the measurable quantities A = 
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for the weighted averaging . 


( 10 ) 


Corresponding errors are calculated using expressions similar to Eq. ([^ for the arithmetic averaging and 
to Eq. 0 for the weighted averaging. Then we used these two ways of averaging to calculate the mean 
values over the hadronic species h, but we do not mix the ways of averaging with each other. 

It is necessary to mention that in some cases before the averaging over the measurable quantities A it was 
necessary to average several sets of the data existing for the same quantity A. Eor example, for y/sj^N = 
4.87 GeV the yields of charged kaons are known for midrapidity and in full 4 vr solid angle. Corresponding 


values (x^/n)y 


were found first for three HG models (HSD, UrQMDl.3 and UrQMD2.1). Then they 


M 

were averaged over two sets of yields (at midrapidity and the full one), and only after these steps they were 
averaged over the types of kaons. This sequence can be found from Table I for the set 1 measured at 
the collision energy = 4.87 GeV. 

Similarly, we performed averaging, if for the same quantity A there were available data in different 
kinematic regions. Eor instance, for the collision energy y/sjsfN = 4.87 GeV the RQMD2.1 model [2T] 
provides the description of distributions of A hyperons for four intervals of longitudinal rapidity in the 
range 2 < y < 3.2 (see Eig. 2 m and A set 1 in Table I). On the other hand, two versions of RQMD2.3 
model were used in |30] to describe the rriT distributions of A hyperons for other four longitudinal rapidity 
in the range 2.2 < y < 3.4 (see Eig. 5 in [30]). Therefore, hrst of all it was necessary to determine the 
QDD of A hyperons over rriT distributions for these models at each rapidity interval yk and then to average 
the obtained values over all rapidity intervals. Then the niT distribution results of A hyperons found for 
two versions of RQMD2.3 model m were further averaged with the QDD for the longitudinal rapidity 
distributions (see Eig. 7 in m and A set 2 in Table I). Such information and the final results can be 
found in Tables I and III-XI for the arithmetic averaging. Eor the weighted averaging Q the subsequence 
of steps was absolutely similar and, hence, for such a way of averaging we give the final results only. In 
the next section we demonstrate a detailed way of finding the QDD for the collision energy y/sNN = 4.87 
GeV. 


III. DETAILS OF AVERAGING PROCEDURE FOR = 4.87 GEV 

To explain in more details the procedure of averaging, in this section we consider it on the example 
of our meta-analysis for the collision energy ^/sJ^ = 4.87 GeV, i.e. for laboratory energy Eiab = 10.7 
GeV. Let us for this purpose analyze the set I from the Table I. To get the QDD of rriT spectra 
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Fig.7, Ref. [SI] 
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Fig.l, 2 Ref. [31] 
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Fig.l, Ref. |37| 
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ARC,RQMD2.1 

Fig. 2 Ref. [21] 

2.75 ± 1.66, 5.74 ±2.1 
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Fig. 4 Ref. [21] 
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Fig. 1 Ref. [31] 
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TABLE I; The QDD provided by HG and QGP models. The 1-st column indicates the particle species, the 2-nd 
one shows results for the quality of rriT spectra description, the 3-rd one shows results for the quality of spectra 
description, while the 4-th one gives results for the QDD of yields at the collision energy y^SNN = 4.87 GeV. In 
some rows there are two values of (x^/n) which correspond to the models and references shown below in the same 
column. For more details see the text. 
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FIG. 1: (Color online) Comparison of transverse mass spectra of K'^ and K mesons at midrapidity from HSD 
(solid curves) and UrQMD 2.0 (dashed curves) for -^Jsmn = 4.87 GeV taken from Fig. 7 of [3l]. The result for 


ix^/n) 




shown in this figure is described in the text. 
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Averaging the above results over types of hadrons one hnds that 


( 12 ) 


= 1.26 ±0.34. At the same 


aaHG 


time the weighted averaging gives us practically the same result {x^In)^^ 


= 1.125 ±0.32. 


waHG 


In the same way we determined the QDD of the longitudinal rapidity distributions of kaons 


(x^/n)^ _ for y^sjsfN = 4.87 GeV using the results of QuarkComb. model 

aaQGP 

set 1 in Table I for details) 


(see Fig. 


and 


ix'^/n) 


Vl 


\K- 


= 4.47 ± 1.22 

QuarkGomb 

= 0.236 ±0.28 

QuarkComb 


ix^/n^ 


K± 

'VL 


{X^ln,y, 

The weighted averaging gives us a different result 
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FIG. 2: (Color online) The longitudinal rapidity spectra of K'^ and K mesons (symbols) reproduced by the 
QuarkComb. model [51]. This is part of Fig. 5 of [M] which was used to determine {x^ln)y^ 
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FIG. 3: (Color online) The yields of and K~ mesons (symbols) reproduced by the following HG models: HSD 
(solid curves), UrQMDl.3 (short dashed curve) and UrQMD2.1 (long dashed curve). The left column corresponds 
to the particle yields measured at midrapidity, while the right column gives the total yields in Itt solid angle. This 


is part of Fig. 2 of m which was used to determine (x^/n)^ 
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From Eqs. (14) and ( |16[ ) we found the corresponding average over two multiplicity sets for positive kaons, 
while Eqs. (15) and allowed us to determine a similar average for negative kaons 
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where in the last step we averaged the results over two kinds of kaons. Repeating the same sequence of 
steps for the weighted averaging of Eqs. (6) and (7), we found 
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about various short hand notations used in Tables I and III-XI, 
which contain only the results of arithmetic averaging. Tables III-XI are given at the end of this work. To 
shorten our remarks inside these tables the sign & is used to demonstrate the fact that the corresponding 
value of (x^/^) is averaged over two sets of data or over the results of models. For example, the notation 
^1^=0^ dvr used in Table I shows that the value of (x^/n) is averaged over the yields measured at 
midrapidity and in the full solid angle. Similarly, the notation ‘HSD & UrQMD2.0’ means that we give the 
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value of (x^/^) averaged over the results of models HSD and UrQMD2.0. In Table I, the notation ‘mr+y’ 
made for A hyperons (set 2) means that we give the value of (x^/n) averaged over rriT- and y-distributions. 
More explaining remarks can be found in the captions of corresponding Tables. 



FIG. 4: (Color online) KTBO plot 1: Comparison of {x^ 

^ aa{A} 


(black symbols and dashed curve) and 


ix^/nY 


a.{h} 


^{A} 


aaHG 


(red symbols and solid curve) as functions of collision energy obtained for the arithmetic averag- 

aaQGP 

ing. The symbols of different hadrons which correspond to the same collision energy are slightly spread around the 


energy value for better perception. The symbols are connected by the lines to guide the eye. 


IV. RESULTS 


Our major results are shown in Figs. and and in Table II. The auxiliary results are given in the 

Tables I, III-XI. As one can see from the KTBO-plot 1, i.e. Fig. there are several matches of the QDDs 

by HG and QGP models. At low collision energies the HG quantity {x ‘^/is smaller than the 

aoHG 

one of QGP models, i.e. QDD is higher for HG generators. Although at energies below 4.87 GeV the QGP 
models are represented by the results of 3FD model, we stress that it would be extremely hard to reach a 
better description of the data by these models compared to the HG ones. It is so, since the latter provide 
almost an excellent description of the analyzed data for 3.1 GeV < ^/snn < 4.2 GeV, as one can see from 
Table H. 
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At the collision energies y/s^N =4.2 GeV =4.87 GeV and y/s^N =12.3 GeV one finds 


ix^/n) 


aa{h} 

aa{A} 


aaHG 


ix^/n) 


aa{h} 

aa{A} 


aaQGP 


(19) 


In the collision energy ranges 4.87 GeV < y/s^N < 12.3 GeV and 12.3 GeV < y/sNN < 17.3 GeV the QGP 
models describe the data essentially better. Therefore, the arithmetic averaging meta-analysis suggests 
that at energies below 4.2 GeV there is hadron phase, while in the region 4.2 GeV < y/s^N < 4.87 GeV 
there is hadron-QGP mixed phase, while at higher energies there exists QGP. Such a picture is well fit into 
the recent findings of the generalized shock adiabat model [151116) . However, the most interesting question 
is how should we interpret the coincidence of two sets of results at the collision energy y/s^N = 12.3 GeV? 




FIG. 5: Schematic pictures of possible locations of the initial states of matter formed in A+A collisions are shown on 
the plane of baryonic density and pressure. Each point on these trajectories (dashed curves) corresponds to a single 
collision energy value. Left panel: As it is argued in the text the possible initial states correspond to the trajectories 
AD or BD as it follows from KTBO-plot 1 for the case of critical endpoint. The trajectory CD is located far from 
the mixed phase region and, hence, it cannot generate the second region in which the QDDs of HG and QGP models 
are equally good. Right panel: In case of the tricritical endpoint the second region in which the QDDs of HG and 
QGP models are equally good may, alternatively, appear due to the second phase transition. 

At first glance it seems that at the collision energy y/s^N = 12.3 GeV the QGP states created by the 
corresponding generators touch the phase boundary with hadron phase. However, one must remember that 
both curves depicted in Fig. have, in fact, finite width defined by the error bars. Taking into account 
an overlap of the curves with finite error bars, one immediately concludes that the overlap region is rather 
wide on collision energy scale, namely it ranges from y/s^N — 10 GeV to y/s^N — 13.5 GeV. Recalling that 
the collision energy width of the mixed phase at low values of ^/snn is below 1 GeV, one may guess that 




























14 


the most probable scenario is that the initial equilibrated states of the matter formed in A+A collisions at 
^ys]s^N — 10 GeV return into the hadron-QGP mixed phase, while at — 13.5 GeV they return from 

the mixed phase into QGP. 

From the present meta-analysis results it is hard to tell, whether the initial equilibrated states of the 
matter formed in A+A collisions at ^/s]^ G [10; 13.5] GeV reach the hadronic phase or they correspond to 
the mixed phase alone. However, such a wide range of collision energy of about 3.8 GeV should correspond 
to something entirely new, compared to the mixed phase located at the range ^/sNN £ [4.2; 4.87] GeV. 
We believe that the energy range y^+viv £ [10; 13.5] GeV may correspond to the vicinity of the critical 
endpoint of the QCD phase diagram. The most probable locations on phase diagram of the initial states 
which formed in A+A collisions are schematically shown in the left panel of Fig. [^in case of the critical 
endpoint. Our main reasons in favor of such a hypothesis are as follows. Recalling argumentation of Ref. 
|18j on the location of (tri)critical endpoint, we have to stress that the vicinity of collision energy 10 GeV 
was independently found in this work as the second entrance into the mixed phase from QGP. On the 
other hand, the lattice QGD results on the baryonic chemical potential fiB dependence of pseudo-critical 


temperature Tps{fJ.B) found from maximum of chiral susceptibility 


or from chiral limit 


show 


that the chemical freeze-out states at ^/sJ^ = 17.3 GeV [HI [12] correspond to the cross-over, but to the 
critical endpoint (see m for an extended discussion). Therefore, we conclude that such a point should 
exist at collision energies below 17.3 GeV, i.e. close to the discussed region of collision energies. 

Of course, alternatively, the second mixed phase region found here at ■,/sj^= 10-13.5 = 11.75 ± 1.75 
GeV may correspond to the second phase transition (chiral), but even in this case the (tri)critical endpoint 
of the QCD phase diagram should be located nearby. Such a situation is depicted in the right panel of 
Fig.[^ 

In order to better determine the energy range of the second entrance inside the mixed phase, let us 
verify our conclusions on the meta-analysis results shown in the KTBO-plot 2 (see Fig.[^. This KTBO-plot 
contains the meta-analysis results for the weighted averaging. Comparing KTBO-plot 1 and KTBO-plot 


waQGP 


^ / 0 / \ wa^h} 

/^La{A} 


waHG 


are essentially 


2 one can see that in many cases the values of 

■’ ' ' walA} 

smaller than the corresponding values obtained by the arithmetic averaging. This is a clear demonstration 
of the fact that the weighted averaging Q and 0 ‘prefers’ the results with minimal error bars, which 


in our case means the smaller value of 


M 


due to Eq. (2). Therefore, one has to be careful with 


this kind of averaging because at each energy there are only a few hadronic characteristics which have low 


value of ln)\ 


M 


Nevertheless, our main conclusions from this plot are qualitatively similar to the case 
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FIG. 6: (Color online) KTBO plot 2: Comparison of 


(black symbols and dashed curve) and 




waHG 


(red symbols and solid curve) as functions of collision energy obtained for the wighted averaging. 

waQGP 

The symbols of different hadrons which correspond to the same collision energy are slightly spread around the energy 
value for better perception. The symbols are connected by the lines to guide the eye. The thin straight lines below 
collision energy 12 GeV demonstrate the intersection region due to finite error bars. 


of the KTBO-plot 1, although there are some quantitative difference. Thus, the equality 


ix^/n) 


wa{h} 

wa{A\ 


waHG 


ix^/n) 


wa{h} 

wa{A\ 


waQGP 


( 20 ) 


holds for the collisions energies y/s^N = 4.87 GeV, y/s^N =10.8-12 GeV and y/sjsrN = 17.3 GeV. 
Accounting for the finite error bars, we found that within the error bars the HG and QGP models describe 
the data equally well in the following region of collision energy: 4.4-5 GeV, 10.8-12 GeV and 14.3-17.3 
GeV. In other words, within the present meta-analysis the region of energy y/sjsfisf = 4.2-4.87 GeV which we 
associated above with the hadron-QGP mixed phase remains almost the same, giving us the overlapping 
region of energies {y/s^N) overlap = 4.4-4.87 GeV. However, instead of a single region of reentering the 
mixed phase now we found two ones. The most peculiar thing is that at the collision energies yT)Tiv = 
11.3-17.3 GeV the weighted averaging of HG models demonstrates somewhat better description of the data 
than the one of QGP models (see Fig. |^. 

How can one understand these results? Our interpretation is as follows. We used two entirely different 
ways of averaging in order to study the stability of the results, and hence the most probable results of 
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Arithmetic averaging 

Weighted averaging 

y/SNN (GeV) 

HG models 

QGP models 

HG models 

QGP models 

3.1 

1.85 ±0.575 

4.1 ±1.84 

0.104 ±0.077 

1.62 ±0.49 

3.6 

0.913 ±0.442 

3.976 ± 1.725 

0.272 ±0.1 

0.94 ±0.395 

4.2 

2.086 ±0.762 

2.616 ± 1.306 

0.1042 ±0.08 

1.14 ±0.46 

4.9 

2.813 ±0.9 

3.57 ±0.608 

1.086 ±0.13 

0.772 ±0.22 

5.4 

6.1 ±0.854 

N/A 

1.733 ±0.346 

N/A 

6.3 

7.24 ±2.06 

1.78 ±0.37 

6.044 ± 1.737 

0.493 ± 0.54 

7.6 

5 ±0.9 

1.08 ±0.27 

2.4 ±0.358 

0.405 ±0.126 

00 

bo 

3.66 ± 1.06 

1.41 ±0.613 

2.355 ±0.208 

0.53 ±0.1 

12.3 

2.29 ±0.9 

2.38 ±0.33 

0.0564 ±0.027 

0.6 ±0.11 

17.3 

10 ± 1.51 

3.475 ± 0.427 

0.38 ±0.0531 

0.434 ±0.12 


TABLE II: The resulting QDD values provided by two ways of averaging. These results of HG and QGP models are 
shown in KTBO-plots 1 and 2. 


KTBO-plots land 2 are the ones which coincide. Therefore, the most probable energy range for the second 
mixed phase region is {^/snn) overlap = 10.8-12 = 11.4 ± 0.6 GeV. 

As an independent check up, let us consider the energy ranges -s/snn = {11.75±1.75; 11.4±0.6; 15.5±1.5} 
GeV of the second mixed phase region which we found from the KTBO-plots 1 and 2 as the results of 
independent meta-measurements. Then applying to them the both ways of averaging, i.e. using Eqs. Q 
and Q, we obtain 

(x/i]^)aa = 12.9 ±0.8 GeV, (= 11.95 ± 0.53 GeV . (21) 


The weighted averaging in (21) gives us nearly the same estimate which we found from the overlapping 
regions in the KTBO-plots 1 and 2, while the arithmetic averaging provides us with a more conservative 
estimate. We suggest to use this value {^SNN)aa = 12.9 ± 0.8 GeV as the most probable estimate for the 
(tri)critical endpoint collision energy. The main reason is that the two other estimates, i.e. s mn) overlap = 
11.4±0.6 GeV and {■s/sNN)wa = 11.95±0.53 GeV, correspond to the second mixed phase region and, hence, 
it is logical to assume that the (tri) critical endpoint energy is located more close to the upper boundary of 
these estimates or even at a slightly higher energy. This is exactly the region provided by the arithmetic 
averaging of the collision energy values. Also such a conclusion is supported by the fact that the endpoint 
was, so far, not found at ^ysNN = 11-5 GeV and ^/s]^ = 12.3 GeV. Recall also the main conclusion of 
RHIG Beam Energy Scan [¥] that the properties of strongly interacting matter created in A±A collisions 
at energies below and above ^ysNN = 11-5 GeV are different. Hence, the value {■s/sNN)aa = 12.9 ± 0.8 
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GeV is the present days best estimate which is provided by the suggested meta-analysis. Hopefully, it can 
be improved further, if one accounts for the RHIC Beam Energy Scan results measured at the collision 
energies ^ysNN = 11-5 GeV and yGjV/v = 19.6 GeV. 

V. CONCLUSIONS AND PERSPECTIVES 


Here we performed the meta-analysis of the QDD of the existing A-|-A event generators without and 
with QGP existence which allow one to extract physical information of principal importance. These 
kinds of generators are, respectively, called the HG and QGP model. A priori we assumed that, despite 
their imperfectness, these models contain the grain of truth on the QGP formation in central nuclear 
collisions which we would like to distillate using the suggested meta-analysis. For each collision energy 


we consider the set of QDD 


of the experimental quantity A of hadron h described by the 


M 


model M G {HG;QGP} as the results of independent meta-measurements. The studied experimental 
quantities include the transverse mass spectra, the longitudinal rapidity distributions and yields measured 
at midrapidity and/or in the full solid angle. In this work we analyzed the strange hadrons only, since 
they provide us with one of the most popular probes both of experimental measurements and of theoretical 


investigations. Using two ways of averaging for the QDD and its error, we were able to extract 

M 

the QDD by two kinds of models at each collision energy. 

Gomparing the results found by these two kinds of models we were able to locate the regions of their 
equal QDD at the collision energies {^/snn) overlap = 4.4-4.87 GeV and {^/snn) overlap = 10.8-12 GeV, 
which we identified with the mixed phase regions. As expected, at center of mass energies below 4.2 GeV 
the HG models ‘work’ better than the QGP ones. As it is seen from the KTBO-plots 1 and 2, in the collision 
energy range = 5-10.8 GeV the HG models fail to reproduce the vast majority of data, while the 

QGP models reproduce data rather well. Therefore, this region we associated with QGP. Unfortunately, 
results of the A-l-A event generators used here did not allow us yet to uniquely interpret our findings at 
the collision energies above 12 GeV. Our educated guess is that the energy range (y^T/vjv)aa = 12.9 ± 0.8 
GeV corresponds to the vicinity of the (tri)critical endpoint. Such a guess is, on the one hand, supported 
by the phenomenological arguments of Ref. |18j on the location of (tri)critical endpoint. On the other 
hand, our hypothesis is also confirmed by the results of lattice QGD |43H46j that the chemical freeze-out 
states at ^/sNN = 17.3 GeV belong to the cross-over region and, hence, the (tri)critical endpoint should 
be located at lower energies of collision. 

Of course, the energies {^ysNN)aa = 12.97 ± 0.92 GeV may correspond to the second phase transition of 
QGD, but even in this case the (tri)critical endpoint of the QGD phase diagram should be located nearby. 
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Perhaps, our findings may help to interpret the conclusions of the RHIC Beam Energy Scan Program [1] 
that below and above ^ysNN = 11-5 GeV one, respectively, probes the different phases of QGP which 
have different properties. We hope that this question could be answered, if the RHIG Beam Energy Scan 
data measured at center of mass energies 11.5 GeV and 19.7 GeV were analyzed by the event generators 
and then reanalyzed using the present approach. We would like to stress that the present work strongly 
supports conclusions of the generalized shock adiabat model nail] that the hadron-QGP mixed phase is 
located between the center of mass energies 4.2 GeV and 4.87 GeV. On the other hand, there is absolutely 
no evidence for the onset of deconfinement at the center of mass collision energy 7.6 GeV, as it is claimed 
in Therefore, we believe that without establishing a much more convincing relation between the 

irregularities discussed in [10] and the deconfinement phase transition one simply cannot consider them to 
be signals of deconfining transition. 

The current set of analyzed models does not cover the full spectrum of existing ones, but already this 
sample was sufficient to obtain interesting results and show the power of the suggested meta-analysis. Also 
it is worth to note that to essentially change the results of present meta-analysis it would be necessary to 
greatly improve the description of existing data either by HG models in the center of mass collision energy 
range 5.4-12.3 GeV or by QGP models at energies below 4.2 GeV. Otherwise, it would be hard to change 
the obtained results, since they are based on the analysis of more than 100 different experimental data 
sets. 

Due to absence of necessary information on the QDD of A-|-A collisions, we were forced to make the 
tantalizing efforts in order to scan both the experimental data and their description by the event generators 
of such collisions from the published hgures. Such efforts could be avoided, if all existing experimental 
data on A-l-A collisions were collected together like it is done in the Particle Data Group for collisions 
of elementary particles. Also the success of present meta-analysis suggests that the results of theoretical 
works which describe data should be presented not only in hgures, but also they should be available as the 
numbers in the form of QDD /n)\ and as the total mean deviation squared per degree of freedom (if 


M 

such a number can be found). Then the hrst of them can be used for the meta-analysis we performed here, 
while the second ones will help us to determine the most successful A-|-A event generators of available data. 
Such an information can be used for planning the collision energy of ongoing experiments and for further 
improvement of the existing A-|-A event generators. We believe that only along this way the heavy ion 
collisions community will be able to solve its ultimate tasks, i.e. to locate the hadron-QGP mixed phase 
and to discover the QCD phase diagram endpoint. Therefore, as the first step in this direction in Tables 
I and III-XI we present the values of QDD used in our meta-analysis for the arithmetic averaging along 
with the final values for the both ways of averaging used in the KTBO-plots 1 and 2 (see Table II). The 
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values of QDD obtained for the weighted averaging will be published elsewhere. 

As the first practical steps we suggest to the RHIC Beam Energy Scan Program to make a few 
measurements in the center of mass collision energy ranges from 3.9 GeV to 5.15 GeV and from 10.8 GeV 
to 15.5 GeV. The hrst of them may lead to a discovery of the hadron-QGP mixed phase, while the second 
one may lead to a location of the QCD phase diagram endpoint. 
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VI. APPENDIX: DESCRIPTION OF THE MODELS COMPARED 

Here we briefly discuss characteristics of the models used in the present meta-analysis. We have to 
apologize for giving a few basic references only. The main criterion for choosing these models was the 
wide acceptance of their results by heavy ion community and a requirement that their results are well 
documented. In our analysis we on purpose included not only cascades, hadronic or partonic, but also 
statistical (SHM) and hydrodynamical (3FD) models, since we believe that the meta-analysis works better, 
if the models of different kinds are included. We consider a model to be of a QGP type, if in it, explicitly 
or implicitly, there appears a dense medium of quarks and gluons (partons), whose properties are entirely 
different from the ones of hadronic medium employed in this model. 

ARC (A Relativistic Cascade) |22[ 08] is a purely hadronic cascade, developed in the Brookhaven 
National Laboratory in the early 90s. ARC is designed for and limited to AGS energies. It does not include 
any potentials, mean field effects, off-shellness, strings or in-medium effects. Hadrons participate in 2 —?■ 2 
elastic and inelastic collisions, resonance formations and decays. Cross-sections and their parameterizations 
are taken from CERN-HERA and LBL compilations as well as from papers [l9]. Collisions are triggered 
using a geometrical criterion: if distance between particles rij < {aij, where aij is corresponding 
total cross-section, then collision happens. Formation time of 1 fm/c is introduced. Apparently, this is a 
HG model. 

RQMD (Relativistic Quantum Molecular Dynamics) [SO] |5T] is a hadron cascade, developed since 1989 
in the Institute of Theoretical Physics of the University of Frankfurt upon Main. It employs nuclear poten- 
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tials, implemented in a Lorentz-covariant manner, at low energies and fills cross-section with interacting 
strings at higher energies. Interactions of energy up to 2-3 GeV mainly proceed through hadronic res¬ 
onances, such as /?, T], K*, (f>, A, N', parametrized by many-channels Breit-Wigner formula. Unknown 
cross-sections are calculated using additive quark model, namely, atot ~ (1 — 0.4xi)(l — 0.4x2), where 
^ 1,2 = jy ’ strange to non-strange quarks ratio in the corresponding hadron. Like in ARC, there is a 
formation time of 1 fm/c and geometrical collision criterion. Hadronic cross-sections are taken from differ¬ 
ent sources, compared to ARC [52]. At higher energies interactions go via strings using phenomenological 
Lund model, extended by possibility for strings to interact, forming so-called ropes (see 1S3| for review). 
RQMD is often believed to be valid at a broad collision energies range from ~50 MeV to hundreds GeV 
per nucleon. 

UrQMD (Ultra-relativistic Quantum Molecular Dynamics) is a successor of RQMD, developed in 
Frankfurt ITP in the late 90s |26j . Most recent version can be run as a hybrid model, but we consider 
only earlier pure hadronic cascade versions. Unlike RQMD, UrQMD employs non-relativistic nuclear 
potentials, that are switched off at relative momentum of colliding hadrons above 2 GeV. Resonance 
table is significantly improved compared to RQMD. Cross-sections are mainly taken from Particle Data 
Group (PDG) and CERN-HERA compilations and are parametrized using multi-channel Breit-Wigner 
ansatz. Experimentally unknown channels, such as resonance-resonance or hyperon-resonance scattering, 
are parametrized using phenomenological additive quark model. As RQMD and HSD, UrQMD employs 
Lund string model at center of mass collision energies y/sArjv > 3 GeV. Unlike RQMD, UrQMD strings do 
not interact. 

HSD (Hadron String Dynamics) is a hadron cascade developed at Giessen in the mid-90s |24l I25j . It 
employs potentials, implemented as a relativistic mean-field. HSD particles propagate off-shell and their 
self-energies are modified in the medium. In the hadronic scattering part HSD is the successor of BUU 
code, developed by the Giessen group [5l|. Collisions happen according to the geometrical criterion, as 
in ARC and RQMD. Cross-sections are taken from similar sources to ARC plus CERN ISR compilation. 
In some cases strange particle cross-sections are re-parametrized to change cross-section behavior in the 
energy regions, where data are not available m- At the center of mass collision energies of nucleons higher 
than 2.6 GeV interaction proceeds via strings |24j, operating according to Lund string model. Strings 
cannot interact: ropes and string fusion like in RQMD are not implemented. Parameters of Lund model 
are changed from energy to energy to match experimental data [25] . 

AGSHIJET |28| is a hadron cascade inheriting from HIJET |56|, which, in turn, comes from a ISAJET 
generator for p + p collisions 1SZ|. All these models were developed at Brookhaven National Laboratory in 
late 80s - early 90s. AGSHIJET operates with strings and uses not the Lund model, but a Field-Eeynmann 
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model, which has a different parametrization for fragmentation functions. HIJET simulates A+A collision 
as a sequence of independent NN collisions, where individual NN collisions are treated using IS A JET. 
Eurther HIJET applies a hadron cascade, similar to ARC, to the obtained products. In the AGSHIJET 
resonance treatment was improved, N* added and Field-Eeynmann model parameters were adjusted to 
reproduce p + p reaction cross-section at laboratory frame momentum below 50 GeV. 

SHM (Statistical Hadronization Model) is not a microscopic, but a pure statistical model of hadron 
gas |271 [58] , which allows one to obtain hadron yields (but not spectra) under the assumption of a sharp 
chemical freeze-out. In this model it is assumed that before chemical freeze-out matter is in thermal 
and chemical equilibrium and after it inelastic reactions are frozen and only resonance decays and elastic 
reactions occur. Thermal yields of hadrons are calculated from the equilibrium ideal gas partition function 
and read as N = f ■ Temperature T, volume V, baryon chemical potential fi and 

strangeness suppression factor 7 ^ are free parameters, adjusted at each collision energy to provide the 
best fit of experimental hadronic yields. Interaction between hadrons is taken into account via adding 
resonances as independent components of ideal gas. 

QuarkComb (Quark Gombination model) [3l| is a phenomenological gluonless model, which generates 
effective quarks with rapidity distribution of = B ^ where B is normal¬ 

ization constant, while a and a are adjustable parameters. Number and flavour of quarks are generated 
according to a phenomenological model of hadronization introduced by Xie and Liu in 1988 [59| and tuned 
to describe hadron production in -|- e“ collisions. Effective (anti)quarks are sorted in rapidity and 
neighbors in the array are combined into mesons or baryons: qqq and its permutations give a meson and 
one quark remains in the array, while qqq gives a baryon. To obtain final spectra hadron resonances are 
decayed. According to our criterion this is QGP model. 

3FD (3-fluid dynamics) is a hydrodynamical model of heavy ion collision, developed by Ivanov and 
Russkikh [35]. It considers ion collision as three fluids: two baryon-rich fluids corresponding to projectile 
and target and the third one, corresponding to produced particles. A following system of coupled equations 
of perfect fluid dynamics 

+ (22) 

= -F^ + F^,, (23) 

+ Fr - f ;, , (24) 


is solved for the equation of state with the first order phase transition between hadronic matter and QGP 


IMIETI. The left hand side of system ( |22[ ) is standard for relativistic hydrodynamics, while its right hand 
side represents interaction between fluids. Since + + T^^) = 0, total energy and momentum are 
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conserved. Freeze-out is performed using a modified Cooper-Frye procedure and, hence, it experiences the 
typical problems of other hydrodynamic models discussed in |60j . After freeze-out the obtained resonances 
are decayed. For considered two phase equation of state we assign 3FD approach to QGP models. 

PHSD (Parton-Hadron-String Dynamics) |38j is the most advanced cascade model based on HSD. 
It combines the HSD hadronic sector with a new model for partonic transport and hadronization, so- 
called DQPM (Dynamical Quasi-Particle Model). DQPM describes QCD properties in terms of single¬ 
particle Green’s functions (in the sense of a two-particle irreducible approximation). Elastic and inelastic 
q + q'^q + q, g + g^g + g^g + g'^g,q + q'^g reaction are included, fulfilling the detailed balance 
condition. Hadronization is occurring continuously based on Lorentz-invariant transition rates. Obtained 
hadron is embedded to the HSD cascade. If the hadron mass is too large, then it is considered as a HSD 
string. Apparently, this is QGP model. 

Core-Corona Model [ID] is a phenomenological approach suggested in 2005 to parameterize and to 
represent experimental data in a simple, but effective way. It was further developed in mi- This model 
is based on the assumption that the fireball created in nuclear collisions is composed of a core, which has 
the same properties as a very central collision system, and a corona, which has different properties and is 
considered as a superposition of independent nucleon-nucleon interactions. The sizes of core and corona 
are defined by the condition that the corona is formed by those nucleon-nucleon collisions in which both 
participants interact only once during the entire process of A-|-A collision. The fraction of single scatterings 
is calculated using straight line geometry as described in the Glauber model |6T]. Since nowadays there 
are no doubts that at sufficiently high energies in central A-|-A collisions the QGP is formed, therefore at 
such energies we regard a core state of this model as an implicit parameterization of QGP state. 
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y/sNN = 3.1 GeV 


TOT-distribution 

rapidity distribution 

Yields 

ix^/n) = 
A(S0) 

N/A 

N/A 

0.856 ± 0.755 4^) 

BSD & UrQMD1.3(2.1) 

Figs. 1, Ref. [3T] 

(y2/n) = 

Ar+ setl 

1.575 ±0.724 

3FD 

Fig. 1, Ref. [33 

N/A 

5.9 ±4.85 

3FD 

Fig. 9, Ref. [36] 

ix^/n) = 

K~ setl 

1.49 ± 0.676 

3FD 

Fig. 1, Ref. [33 

N/A 

7.43 ± 5.45 

3FD 

Fig. 9, Ref. [36] 

ix^/n) = 

K+ set 2 

1.35 ±0.476 

BSD & UrQMD2.0 

N/A 

8.246 ± 3.316 ^1 ^ 

dy \y=0 

BSD & UrQMD1.3(2.1) 

Fig. 1 or 2, Ref. |3I| 

(xV«) = 

K~ set 2 

1.7 ±0.51 

BSD & UrQMD2.0 

N/A 

0.0714 ± 0.31 n 

dy \y=0 

BSD & UrQMDl.3(2.1) 

Figs.l or 2, Ref. jS] 

(xV«) = 

1.525 ±0.35 

BSD & UrQMD2.0 

Fig. 7, Ref. [U] 

N/A 

N/A 


TABLE III: Same as in Table but for collision energy ^/SNN =3.1 GeV. 
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^snn = 3.6 GeV 


mr-distribution 

rapidity distribution 

Yields 

(xV«) = 

A(E0) 

N/A 

N/A 

0.862 ± 0.76 4^) 

HSD & UrQMDl.3(2.1) 

Figs. 1, Ref. [31] 

(xV«) = 

K~^ set 1 

5.53 ±2.74 

mT-distr.±Yield:3FD 

mT: Fig. 1, Ref. [37] ;Yield:Fig.9,Ref. [35] 

N/A 

2.13±1.7f/^| J 

dy \ y=0 j 

HSD & UrQMDl.3(2.1) 

Fig. 1 or 2, Ref. [3T] 

ix^/n) = 
A'+ set 2 

1.17 ±0.46 

HSD & UrQMD2.0 

Fig. 7, Ref. [311 

N/A 

N/A 

(xV«) = 

K~ setl 

0.453 ±0.286 

HSD & UrQMD2.0 

Fig. 7, Ref. [311 

N/A 

0.0824 ± 0.332 J 

V dy \y=0j 

HSD & UrQMDl.3(2.1) 

Fig. 1 or 2, Ref. [3T] 

(xVn) = 

K~ set2 

0.494 ± 0.424 

3FD 

Fig. 1, Ref. [37] 

N/A 

4.35 ±4.17 

3FD 

Fig. 9, Ref. [35] 


^snn = 4.2 GeV 


mT-distribution 

rapidity distribution 

Yields 

(xVn) = 
A(E0) 

N/A 

N/A 

3.35 ± 1.47 47r') 

V dy ly=0 J 

HSD & UrQMDl.3(2.1) 

Fig. 1, Ref. [31] 

(xVn) = 

0.81 ±0.3 

HSD & UrQMD2.0 

Fig. 7, Ref. [31] 

N/A 

0.8343 ±0.75 (f 

HSD & UrQMDl.3(2.1) 

Fig. 1 or 2, Ref. [31] 

(xVn) = 

K+ 

3.367 ± 1.22 

3FD 

Fig. 1, Ref. [37] 

N/A 

5.7 ±4.78 

3FD 

Fig. 9, Ref. [ig 

(xV«) = 

K- 

0.73 ±0.54 

3FD 

Fig. 1, Ref. [37] 

N/A 

.665 ± 1.63 

3FD 

Fig. 9, Ref. [3g 


TABLE IV: Same as in Table but for collision energies = 3.6 and 4.2 GeV. 
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y/sNN = 5.4 GeV 


m-r-distribution 

rapidity distribution 

Yields 

ix^/n) = 

N/A 

8.17±2 

BSD 

Fig. 14, Ref. [15] 

N/A 

{x^/n) = 

A set 1 

N/A 

0.45 ±0.6, 1.64 ±1.14 

ARC(SiSi(Pb)) 

Figs. 4-5, Ref. [H] 

1.1 ±2.1, 1.35 ±2.3 

ARC(SiSi(Pb)) 

Ref. [H] 

{x^/n) = 

A set 2 

N/A 

2.24 ±1.3, 6.1 ±2.2 

AGSHIJET-N*(SiSi(Pb)) 

Figs. 4-5, Ref. [H] 

8.9 ±6, 10.2 ±6.4 

AGSHIJET-N*(SiSi(Pb)) 

Ref. [H] 

{x^/n) = 

N/A 

11.7 ±3, 22 ±4 

ARC(SiSi(Pb)) 

Figs. 2-3, Ref. [29] 

1±2, 0.3±1 

ARC(SiSi(Pb)) 

Ref. [29] 

{X^ In) = 

A set 3 

N/A 

1.06 ±0.9, 2.4 ± 1.4 

AGSHIJET-N*(SiSi(Pb)) 

Figs. 2-3, Ref. [29] 

3.5 ±3.7, 7±5 

AGSHIJET-N*(SiSi(Pb)) 

Ref. [29] 


TABLE V: Same as in Tablebut for the collision energy = 5.4 GeV. The notation ‘ARC(SiSi(Pb))’ used 

for A hyperons (set 1) of y-distribution means that at preceding row of this column we give two values of {x^ In) 
which correspond to Si+Si and Si+Pb collisions calculated for the event generator ARC. 
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^snn = 6.3 GeV 


TOr-distribution 

rapidity distribution 

Yields 

ix^/n) = 

1.215 ±0.73 

mT-distr.±Yield:3FD 

mr :Fig. 1 ,Ref. [37] ; Yield:Fig. 9,Ref. [35] 

2.28 ±0.522 

QuarkComb. model 

Fig. 2 Ref. |31| 

N/A 

(y2/n) = 

(t> 

N/A 

1 ± 0.72 

QuarkComb. model 

Fig. 2 Ref. |31| 

2.91 ±3.41, 8.2 ±5.73 

SHM, UrQMD 

Fig. 17 Ref. [31] 

ix^/n) = 

A 

N/A 

5.645 ± 1.5 

QuarkComb. model 

Fig. 2, Ref. [m 

N/A 

ix^/n) = 

A 

N/A 

0.266 ±0.516 

QuarkComb. model 

Fig. 2, Ref. (an 

N/A 

ix^/n) = 

N/A 

0.277 ±0.39 

QuarkComb. model 

Fig. 2, Ref. [51] 

N/A 


TABLE VI: Same as in Table but for the collision energy = 6.3 GeV. 














29 



^snn = 7.6 GeV 


mr-distribution 

rapidity distribution 

Yields 

ix^/n) = 

1.51 ±0.53 

3FD 

Fig. 1, Ref. ISZl 

0.6 ±0.265 

QuarkComb. model 

Fig. 2 Ref. [31] 

N/A 

ix^/n) = 

</> 

N/A 

1.413 ±0.752 

QuarkComb. model 

Fig. 2, Ref. |31| 

0.567 ±1.5, 1.44 ±2.4 

SHM,UrQMD 

Fig. 17, Ref. im 

ix^/n) = 

A 

N/A 

0.4 ±0.4 

QuarkComb. model 

Fig. 2, Ref. im 

N/A 

ix^/n) = 
A 

N/A 

0.363 ±0.5 

QuarkComb. model 

Fig. 2, Ref. |34| 

N/A 

ix^/n) = 
2± 

N/A 

2.16 ±0.856 

QuarkComb. model 

Fig. 2, Ref. |34| 

N/A 


TABLE VII: Same as in Table but for the collision energy = 7.6 GeV. 
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y/sNN = 8.8 GeV 


mr-distribution 

rapidity distribution 

Yields 

{X^ln) = 

set 1 

4.03 ±0.9 

3FD 

Fig.I,Ref. [37j 

2.76 ±0.56, 3.826 ±0.654 

BSD, PHSD 

Fig. 15 Ref. [Sni 

0.1 ±0.3 

Core Corona 

Fig. 10 Ref. [32] 

{X^jn) = 

set 2 

2.33 ±0.7 

PHSD 

Fig. 16 Ref. [Sni 

0.98 ±0.333 

QuarkGomb. model 

Fig. 2 Ref. [Ml 

10.7 ±0.7 

y±Yield: HSD & UrQMD2.3 

Figs. 6-7,10 Ref. [32] 

ix^/n) = 

N/A 

1.8 ±0.846 

QuarkGomb. model 

Fig. 2, Ref. [Ml 

0.94 ± 1.94, 8.614 ±5.9 

SHM,UrQMD 

Fig. 17, Ref. [32] 

(xV«) = 

A(E0) 

N/A 

1.51 ± 0.78, 0.185 ± 0.27; 0.6 ± 0.5 

HSD, PHSD; QuarkGomb. model 

Fig. 17 Ref. |39|; Fig. 2 Ref. [33] 

N/A 

(xV^) = 

A 

N/A 

2.535 ± 1.126, 0.586 ±0.54 

PHSD, QuarkGomb. model 

Fig. 18 Ref. [39], Fig. 2 Ref. [Si] 

N/A 

(xV«) = 

N/A 

1.14 ±0.7, 0.95 ±0.667 

PHSD, QuarkGomb. model 

Figs. 19-20 Ref. [32], Fig. 2 Ref. [M] 

N/A 


TABLE VIII: Same as in Table but for the collision energy y'sjvw = 8.8 GeV. The notation ‘y+Yield’ used for 
iL± mesons (set 2) means that we show the value of {x^/n) averaged over the y-distribution and the particles yield, 
first, and then the results obtained for positive and negative kaons are arithmetically averaged again. Note that the 
results taken from Fig. 10 of Ref. [32] include the most central collisions with N-^iound > 300. 
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^/sNN = XZ.'d GeV 


m-r-distribution 

rapidity distribution 

Yields 

ix^/n) = 

1.26 ±0.4, 5.05 ±0.75 

PHSD, 3FD 

Fig. 16 Ref. [55]. 
Fig. 1, Ref. [37] 

1.3 ±0.4, 2.2 ±0.5; 0.63 ±0.27 

HSD, PHSD; QuarkComb. model 

Fig. 15 Ref. [53]; 

Fig. 2 Ref. [31] 

N/A 

ix^/n) = 

N/A 

6.16 ± 1.326 

QuarkComb.model 

Fig. 2 Ref. [M] 

1.6 ±2.53, 5.264 ±4.6 

SHM,UrQMD 

Fig. 17 Ref. [32] 

ix^/n) = 
A(E0) 

N/A 

0.553 ± 0.526, 0.076 ± 0.2; 1.7 ± 0.923 

HSD,PHSD; QuarkComb. model 

Fig. 17 Ref. [39]; Fig. 2 Ref. [34] 

N/A 

ix^/n) = 
A(S0) 

N/A 

3.2 ± 1.46, 0.24 ±0.4 

PHSD, QuarkComb.model 

Fig. 18 Ref. [33], Fig. 2 Ref. [M] 

N/A 

ix^/n) = 
2± 

N/A 

0.74 ±0.6, 1.453 ±0.7 

PHSD, QuarkComb.model 

Figs. 19-20 Ref. [33], Fig. 2 Ref. [M] 

N/A 


TABLE IX: Same as in Table but for the collision energy s/snn = 12.3 GeV. 
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y/sNN = 17.3 GeV 


TOT-distribution 

rapidity distribution 

Yields 

{x^/n) = 

set 1 

1.83 ±0.2 

5 versions of HSD& 

UrQMD2.0(2.1) 

Figs. 7,8,10,12 Ref. [S] 

25.03 ±0.81 

2 versions of HSD& 

UrQMD2.0(2.1) 

Figs. 9,11 Ref. [31] 

0.6±0.44(f 4^) 

HSD & UrQMD1.3(2.1) 

Fig. 1 Ref. [31] 

(xV«) = 

set 2 

2.675 ±0.667, 0.724 ±0.348 

3FD, PHSD 

Fig.l, Ref. [37]: 

Fig. 16 Ref. [51] 

4.1 ± 0.665, 18.11 ±0.7 

PHSD, HSD & UrQMD2.3 

Fig. 15 Ref. [31], 

Figs. 6-7 Ref. [22] 

1.3 ±0.81, 0.0165 ±0.19 

HSD&UrQMD2.3, 

GoreCorona 

Fig. 10 Ref. [22] 

(xV«) = 

A(E0) set 1 

N/A 

1.53 ±0.78, 1.314 ±0.725 

HSD,PHSD 

Fig. 17 Ref. [31] 

0.66 ± 0.66 4^) 

HSD & UrQMD1.3(2.1) 

Fig. 1 Ref. [31] 

ix^/n) = 
A(E0) set 2 

N/A 

N/A 

0.512 ± 1.43, 0.44 ± 1.327 

HSD, PHSD 

Fig. 21 Ref. [31[ 

ix^/n) = 
A(S0) 

N/A 

10.87 ± 2.33, 18.57 ±3.05 

HSD, PHSD 

Fig. 18 Ref. [31] 

1.4 ± 2.37, 0.0415 ±0.4 

HSD, PHSD 

Fig. 21 Ref. [31] 


TABLE X: Same as in Table but for collision energy ■s/sn'n = 17.3 GeV. The notation ‘5 versions of 
HSD&UrQMD2.0(2.1)’ shown for meson (set 1) means we give the values of (x^/n) averaged over 5 types 
of HSD parameterization (see Ref. [21]) and over 2 versions of UrQMD (see ibid.). Note that the results taken from 
Fig. 21 of Ref. [31] and from Fig. 10 of Ref. [22] include the most central collisions with N^ound > 300. 
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y/sNN = 17.3 GeV 


mT-distribution 

rapidity distribution 

Yields 

{x^/n) = 

0 

N/A 

N/A 

44.6 ± 9.44 

SHM & UrQMD 

Fig. 17 Ref. |35| 

{x^/n) = 

N/A 

9.19 ± 1.75, 1.98 ±0.812 

HSD,PHSD 

Figs. 19-20 Ref. [31] 

0±0; 5.513 ± 3.32, 2.02 ± 2.01 

UrQMD with reduced masses; 

HSD,PHSD 

Fig. 2 Ref. [35]; Fig. 22 Ref. (33] 

{x^/n) = 

A +A 

N/A 

N/A 

0.34 ± 1.166 

UrQMD with reduced masses 

Fig. 2 Ref. [32] 

{x^/n) = 

N/A 

N/A 

3.65 ±3.82 

UrQMD with reduced masses 

Fig. 2 Ref. [35] 


TABLE XI: Continuation of Table for collision energy = 17-3 GeV. Note that the results taken from Fig. 

22 of Ref. [3H] and from Fig. 2 of Ref. [35] include the most central collisions with N^ound > 300. 













